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Abstract. The Falicov-Kimball model is a simple quantum lattice model that describes 
light and heavy electrons interacting with an on-site repulsion; alternatively, it is a 
model of itinerant electrons and fixed nuclei. It can be seen as a simplification of the 
Hubbard model; by neglecting the kinetic (hopping) energy of the spin up particles, 
one gets the Falicov-Kimball model. 

We show that away from half-filling, i.e. if the sum of the densities of both kinds 
of particles differs from 1, the particles segregate at zero temperature and for large 
enough repulsion. In the language of the Hubbard model, this means creating two 
regions with a positive and a negative magnetization. 

Our key mathematical results are lower and upper bounds for the sum of the lowest 
eigenvalues of the discrete Laplace operator in an arbitrary domain, with Dirichlet 
boundary conditions. The lower bound consists of a bulk term, independent of the 
shape of the domain, and of a term proportional to the boundary. Therefore, one lowers 
the kinetic energy of the itinerant particles by choosing a domain with a small boundary. 
For the Falicov- Kimball model, this corresponds to having a single 'compact' domain 
that has no heavy particles. 

Keywords: Falicov-Kimball model, segregation, phase separation. 
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1. Introduction 
1.1. The Falicov-Kimball model. Introduced thirty years ago to describe the semicon- 



ductor-metal transition in SniBg and related materials |FK|, the Falicov-Kimball model 
is a simple lattice model with rich and interesting properties. The system consists of two 
species of spinless electrons with different effective masses: one species has infinite mass (so 
the particles do not move — we call them 'classical particles'), while the second species 
represents itinerant spinless electrons whose kinetic energy is represented by a hopping 
matrix. The Hamiltonian in a finite domain Jl C Z"^ is 



HKiiw.}) 



Yl 4:Cy + 2dY,na: + Uj2 



WxTl 



(1.1) 



x,y£fi 
\x-y\=l 



xGfl 



x£fl 



J 



n. 



CxCx'j 



Here, Cx, Cx, denote creation, annihilation operators of an electron at site x; 
Wx = 0, 1 is the number of classical particles ('heavy electrons') at x, and U ^ is an 
on-site repulsion between the two species of particles. Hq{{wx}) represents the energy of 
the electrons under a potential Uwx- The term 2d^nx in ( |1.1[) is for convenience only. 
It makes H^ positive, and this term only adds 2d times the electron number, N. At zero 
temperature, one is typically interested in the configurations of classical particles that 
minimize the ground state energy of the electrons. 



The model was reinvented in [KL| as a simplification of the Hubbard model, by neglect- 
ing the hoppings of electrons of spin | , say. This simplification changes the nature of the 
model somewhat, mainly because the continuous SU(2) symmetry is lost. Connections 
between the two models are therefore not immediate; however, the greater knowledge 
obtained for the Falicov-Kimball model may help in understanding the Hubbard model. 

Rigorous results in | KL ] include a proof that equilibrium states display long-range order 
of the chessboard type when both species of particles have density 1/2; this holds for all 
dimensions greater than 1 and for allU ^ (including U < 0), provided the temperature 
is low enough. The model is refiection positive under a suitable magnetic field, or when 
the electrons are replaced by hard-core bosons; this property can be used to establish 
long-range order [ |MP| ]. Perturbative methods allow for an extension of these results for 
large U and small temperature, see p^, |MM , DFF|. Absence of long-range order when 
the inverse temperature /3 is small, or PU is small, was also established in |KL|. 

One may increase the density of one species and decrease the density of the other species 
while maintaining the half- filling condition, namely that the total density is 1. (However, 
as was shown in [KL], the lowest energy is achieved when both species have density 1/2.) 
The one-dimensional case was considered in | Lem | ; if classical particles and electrons have 
respective densities ^ and 1 — ?, the ground state is the 'most homogeneous configuration' 
for U large enough; this configuration is periodic with a period no greater than q. Away 
from half-filling the particles segregate: classical particles occupy one side of the chain, 
leaving room for electrons on the other side. There are several results in 2D. Gruber et. 
al. | GJL[| performed a 1/U expansion and found periodic phases when the density of the 
classical particles is 1/3, 1/4, 1/5. This was made rigorous by Kennedy [Ken|. These 
results are reviewed in [PM |. The knowledge of the 2D phase diagram for large U was 
further extended in | Ken2| , Ha^ , |HK| ]. New ground states for various rational densities 
were uncovered; for some densities they are periodic, but there are also densities where 
coexistence of configurations with different periods have minimum energy. The results are 
summarized in Fig. 1 in |HK|. Finally, it was understood in [DMN] that the 111 interface 
is stable, due to the effective interactions between the particles. 
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1.2. Away from half-filling. The purpose of our study is to explore the situation away 
from half-filling; i.e., we take the total density to differ from 1. For any density away from 
half- filling we prove that the ground state is segregated for U large enough. When U = oo, 
the ground state is segregated for all densities (at half filling all configurations have same 
energy, including segregated and periodic ones). Hole-particle symmetries for both species 
of particles ]KL[ imply that the results for positive U and densities (ne,nc) of electrons 
and heavy (classical) particles, transpose to (a) positive U and densities (1 — Ue, 1 — ric), 
and (b) negative U and densities (ne, 1 — Uc) or (1 — n^,nc). For simplicity, we take the 
total density rig -|- n^ to be strictly less than 1. 

We start our study by taking the limit [/ — > cx). Electrons are described by wave 
functions that vanish on sites occupied by the classical particles, and the question is to 
find the arrangement of classical particles that minimizes the energy of the electrons. This 
amounts to minimizing the sum of the lowest eigenvalues of the discrete Laplace operator 
with Dirichlet boundary conditions. This is explained in Section |2|, where it is shown that 
the energy per site of A^ electrons in a finite domain A C Z with volume |A|, is bounded 
below by the energy per site of the electrons in the infinite lattice with density n = A^/|A|. 

One can refine this lower bound by including a term proportional to |5A|, the volume 
of the boundary dh. of A (Section y). This implies that the configuration of the heavy, 
fixed electrons that minimizes the ground state energy of the movable electrons has, more 
or less, one large hole with relatively small perimeter. Thus the movable particles are 



separated from the fixed ones. This behavior was conjectured in [FF| and is opposite to 



the checkerboard configuration, in which both kinds of particles are inextricably mixed. 



Segregation was shown to occur in the ground state of the d = 1 model in | Lem | , and of 



the d = oo model in |FGM|. The present paper proves that this holds for all dimensions, 
and in particular for the relevant physical situations d = 2 and d = 2>. 

Segregation is more difficult to understand on a heuristic level than the chessboard 
phase. The latter is a local phenomenon that results from effective interactions between 
nearest neighbor sites, while the former is a global phenomenon involving extended wave 
functions. This remark should also apply to the Hubbard model, for which antiferromag- 
netism is much better understood than ferromagnetism. 

The fact that the sum of the lowest A^ eigenvalues of the Laplacian in a domain of volume 
|A| is bounded below by the infinite volume value at the same density is not unexpected 



and holds also in the continuum. Indeed, the original idea, due to Li and Yau |LY] (see also 



[ LL | Sections 12.3 and 12.11), was demonstrated in the continuum, and we only adapted 
it to the lattice context. However, the fact that the error term is proportional to |9A|, the 
area of the boundary, is a completely different story. Its proof, at least the one given here, 
is complicated. More to the point, such a bound does not hold in the continuum. One can 
easily construct a continuum domain with finite volume |A|, but with \dA\ = oo, and for 
which all eigenvalues are finite. 

Taking U large instead of infinite decreases the energy, but the gain is at most pro- 
portional to |9A|, as explained in Section |^. Therefore the infinite-C/ segregation effect 
outweighs this gain and the particles are still separated (at least when the heavy particle 
density is far enough from 1/2). Finally, an upper bound derived in Section ^ shows that 
the energy of electrons in A really consists of a bulk term independent of the shape of A, 
plus a term of the order of the boundary. 



These results are summarized in Theorem 1.1 below. One can use them to discuss 



the electronic free energy at inverse temperature /3, for a fixed configuration of classical 



particles, see Theorem 1.2, The conclusion of this paper involves a discussion of first-order 
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phase transitions at finite temperature, of wliat liappens wlien classical particles have a 
small hopping term, and of the possible links with ferromagnetism in systems of interacting 
electrons with spins. 

In order to present the main result of this paper, we need a few definitions. For k G 

(— vr, vr]'^, we set 

d 
efc = 2(i-2^cos/ci. (1.2) 



i=l 



The energy per site e(n) of a density n of free electrons in the infinite volume Z is 

1 



where the Fermi level ep = ^pin) is defined by the equation 

We can specify the configuration {wx)xez<i o^ classical particles by the domain A C Z'^ 
consisting of those sites without particles (holes), that is, w{x) = 1 if x ^ A and w{x) = 
if X G A. Let hj[ denote the one-particle Hamiltonian whose action on a square summable, 
complex function tp on Z is 

[h'iip]ix) = - Yl V>{y) + ^dipix) + UXA4xMx). (1.5) 

y,\y-x\=l 

Here, XA':(a^) is the characteristic function that is 1 if x belongs to the complement A'^ of 
A, and is if x ^ A'^. We define Ej[j^ to be the ground state energy of N electrons for 
the configurations defined by A, i.e. 

TV 
El^= inf V((p.,/i^(^.), (1-6) 

{ipi,...,ipN}~^ 

where the infimum is taken over N orthonormal functions, i.e. {ipi^ipj) = 6ij. There exist 
normalized minimizers if the Fermi level is below U; they are not identically zero inside 
A, and decay exponentially outside. 

Notice that E'^j^ is increasing in U, since {ip, hJ^ip) is increasing in U for any p. 

We define the boundary by dh = {x G A : dist (x, A'^) = 1}, where A'^ is the complement 
of A, i.e., the points in 1/ not in A. The following theorem summarizes the results obtained 
in this paper. It contains upper and lower bounds for the ground state energy. We set 
n = iV/|A|. 

Theorem 1.1. There are functions a{n) > and ^{U) with limjy^oo U^{U) = 8d'^ , such 
that for all finite domains A, 

{2dn - e{n))\dA\ ^ eIj^ - \A\e{n) ^ {a{n) - j{U))\dA\. 
Furthermore, a{n) = a(l — n), and for n ^ \Sd\/{'^T^Y , it can he chosen as 

2'^-3 , , 2 



'^(^) — . rJTT;^ 



1+ 



d . 



Here, \Sci\ is the volume of the unit sphere in d dimensions. 
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An explicit expression for "f{U) can be found in Proposition |4.1| . The lower bound 
a{n) vanishes when n = (no itinerant electrons) or n = 1 (fully occupied lower band). 
Theorem |l.l| is relevant only when a{n) > 7(C/), that is, sufficiently away from half-filling 
(depending on U). The theorem states that the 'good' configurations A for which electrons 
have low energy must have small boundaries. As a consequence, the system displays phase 
separation in the ground state. 

The upper bound is symmetric under the transformation n i— > 1 — n due to a symmetry 
of the Hamiltonian, and it is saturated for f7 = oo by configurations with isolated holes. 
Indeed, in this case the eigenstates consist of 5 functions on the holes, with eigenvalues 
equal to 2d, and dA = A. 

The lower bound is first explained in Section g for U = oo and without the term 



involving the boundary. The latter requires more effort and is derived using Lemmas ^^ 



is proved in Section |5|. 



3.5| in Section y. Proposition |4.l| then extends it to the case of finite U. The upper bound 



Theorem LI is described in | FLU |, which also reviews the rigorous results obtained so 



far for the Falicov-Kimball model. 

1.3. Electrons at low temperature. It is natural to consider the situation at pos- 
itive temperature. The relevant object is the Gibbs state obtained by averaging over 
the configurations of classical particles, and by taking the trace of the Gibbs operator 
exp{— /Jf^Q ({wa;})}. We expect the system to display a first-order phase transition in the 
grand-canonical ensemble; densities of both types of particles should have discontinuities 
as functions of the chemical potentials. But a rigorous treatment of this phase transition 
is beyond reach at present. However, we do obtain some properties of the system when 
the configuration of the classical particles is fixed, and the electrons are at positive tem- 
perature. Namely, one can extend the estimates of the ground state energy to estimates 
of the electronic free energy. The results are described in this section, and their derivation 
can be found in Section |^. 

Let us consider a box fi with periodic boundary conditions. The configuration of classi- 



cal particles is specified by the set of holes A C il (later, in Corollary l.S, we shall average 
over A). With fi being the chemical potential, the grand-canonical electronic free energy 
(equal to — |il|//? times the pressure) is 

i"n,A(/3,/«) = -^logTV exp{-mK,A + Pf^Nn}. (1.7) 



Here, H^ ^ = H^{{wx}) as defined in (LI), Nq = J2x&n''^x is the number of electrons in 



Q, and the trace is in the Fock space of antisymmetric wave functions on $7. 

A simple 'guess' for Fq ^ is obtained by considering independent electrons, which are 
either in A or else in A*^. In the latter case the effective chemical potential is // — C/. Our 
'guess' would then be 

F}{^^{P,fi) « |A|/(/3,/z) + m - |A|)/(/?,/x - U), (1.8) 

where /(/?, /i) is the free energy per site for free electrons: 

^(^'^) = -^7^ / dA;log(l+e-^(^'=-'^)). (1.9) 



Formula ( |1.8D is, indeed, correct when U is large — in the sense that the error is 
proportional only to |i9A|. More precisely, 
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Theorem 1.2. There are functions a{(3,fj,) > with lim^^oo q(/3, /i) > if < fi < Ad, 

and 7(f/) with Yvmu^ao U^{U) = \Qd^ + 2 "^^d^, such that for all finite domains Q. and 

Cd,,\dK\+C',^^M''--^ ^ Fj,%(/3,/i)-{|A|/(/3,/i) + (|0|-|A|)/(/3,/.-t/)} 

^ (a(/3,/x)-7(C/))|M|, 



with 



AirVd „ , ,\ 1 



^-^.^=(^^ + 242^ + 1 



AirVd 1 



The term \Q\^ ^ on the left side is not exactly proportional to |5A|. However, we have 
in mind that |A| and \^\ are comparable, in which case |0| ^a is no greater than \dA\ (up 
to a factor). 

Notice that the upper bound vanishes as // — > — oo, i.e. when the density tends to 0. In 
the limit C/ — > oo. Theorem |1.2| takes a simpler form, namely 

C^JdA\ ^ FEj°°iP,ti)-\A\f{(3,i^) ^ a{P,i^)\dA\. (1.10) 



This extension of Theorem 1.1 to the case of positive (electronic) temperatures is ex- 
plained in Section ^. The lower bound follows from Propositions |6.1| and 3.2, while the 



upper bound is stated in Proposition 6.3. 

Our next step is to find upper and lower bounds for the total grand-canonical 'free 
energy' by averaging over A (i.e., averaging over the positions of the classical particles). 
This can be done with the aid of the Ising model free energy fisingiP, h), 

/ismg(/5,/i) = -^ 1™ Tj^ Yl expj-/? Y^ s^Sy-phY^xj, (1.11) 

\x-y\=l 

where the sum is over configurations of classical spins on 17. 

Corollary 1.3. If U is large enough (so that a{p,fj,) — "){U) > Q), we have Ising hounds 
for the full free energy, 

\[f{P,IJ.) + f{l3,li-U)] + \a + fJi^^^{i-^afi,^[f{l3,,,)-f{l3,,,-U)]) 

^ -2 [/(/?, fi)+f{f3,fi- U)] + idQ,^ + ^Cd,^,fMnSCd,>.l3, c^ [/(/3, /i) - /(/3, /x - [/)]) 
where a = a{(3,fj,) — ^{U). 

The proof can be found at the end of Section y. 

Another consequence of Theorem ^^ concerns the equilibrium state; namely, it allows 
for a precise meaning of segregation. We consider the probability that sites x and y are 
both occupied by classical particles, or both are unoccupied. Namely, we consider 

,. > _ EAcf7;^.=^,exp{-/3F^,^(/3,^)} 

EAcQexp{-/3F//^^(^,^)} 
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where the sums are over subsets A of il such that |A| = [(1 — nc)|ri|] {[z] denotes the 
integer part of z G M). The restriction Wx = Wy means that either both x and y belong to 
A, or both belong to the complement of A. 

Segregation means that up to a small fraction of sites that are close to the boundary 
between classical particles and empty sites, any two sites at finite distance are either both 
hosts of a classical particle, or are both empty. The fraction of sites close to the boundary 
vanishes in the thermodynamic limit. Hence we expect that 

lim lim lim (<5^„,,^Jq = 1, (1.13) 

p— >oo \x—y\~^oo |n|— ►oo 



but we are unable to prove it. Notice that using Theorem 1.1, one can conclude that 



\x—y\^oo\\l\—*oo p— ►oo 

Indeed, taking the limit of zero temperature at finite volume, the sum over A becomes 
restricted to the ground state configuration(s), whose boundary fraction |5A|/|A| tends to 
zero in the thermodynamic limit. 



We can however take advantage of Theorem 1.2 to obtain a result that is better than 



(1.14): 



Corollary 1.4. IfU is large enough (depending on jjl and d only), the ground state of the 
Falicov- Kimball model displays segregation, in the sense that 

lim lim lim ((5^,,t„^)Q = 1. 

|a;— y|— >oo p— >oo |SZ|— +oo 

The proof of this corollary can be found at the end of Section g. 

2. The discrete Laplace operator in a finite domain 

2.1. Basic properties. We start our investigations by taking the limit U -^ oo. Let us 
denote h\ = /i^^°°, the corresponding Hamiltonian, which acts on functions ip E -^^(A) 
as follows: if x G A, 

[hA^]ix) = - Yl Viy) + '^d^ix). (2.1) 

y&A,\x-y\=l 

Some observations can readily be made, that will be useful in the sequel. For cp £ L^(A), 
one has the following formula, 

{(p,hA(p)= Y^ \(p{x) - (p{y)\ , (2.2) 

{x,y}:\x-y\=l 

where the sum is over all x,y £ Z'^, with the understanding that ^{x) = for x ^ A. 
Equation ( |2.2| ) takes a simple form because of the diagonal term in h\ . The effect of the 
Dirichlet boundary conditions appears through a term YlxedA l¥'(^)Pi ^^^^ is due to pairs 
{x, y} with X £ A and y ^ A. 

The matrix h\ is self-adjoint, and ( ^.2| ) shows that h\ ^ 0. Its spectrum has a 
symmetry. Let p be an eigenvector with eigenvalue e, and (p be the function pix) = 
{—iy^'ip{x); one easily checks that (p is also an eigenvector, with eigenvalue {4d — e). 
The spectrum is therefore contained in the interval [0, 4(i], and is symmetric around 2d. 
Furthermore, one has 

EA,\A\-N = 2d\A\il-2n) + EA^N, 

e(l-n) = 2(i(l-2n) + e(n). (2.3) 
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This allows to restrict ourselves to the case n ^ 2 j indeed, existence of a lower bound 
for a density n implies a lower bound for the density 1 — n. This symmetry holds only for 

U = oo. 

2.2. The bulk term. We are looking for a lower bound for the sum E^^n of the first 



N eigenvalues of /ia- This problem was considered by Li and Yau |LY] for the Laplace 
operator in the continuum. Let A C M"^ be a bounded domain. They prove that the sum 
Sn of the first A^ eigenvalues of the Laplace operator with Dirichlet boundary conditions 
is bounded below, 

Sn > {27rfj^\Sd\--^N'+i\A\~i, (2.4) 

where 15*^1 and |A| are the volumes of respectively the d-dimensional sphere and of A. 



The corresponding inequality in the discrete case — our Theorem 1.1 without the bound 



ary correction — constitutes the heart of this paper, and we explain below the proof of Li 



and Yau; see also |LL|, Theorem 12.3. 

It is useful to introduce the Fourier transform for functions in L^(Z'^), that is, on the 
whole lattice. A function tp G -^^(A) can be considered an element of L^(Z'^) by setting 
(p{x) = outside A. This Fourier transform will lead to the electronic energy density for 
the infinite lattice, which is the bulk term for the energy of electrons in A. 

The Fourier transform of a function 93 is defined by 

(^(fc)= ^(^(x)ei'=^ A;e[-^,7r]^ (2.5) 

and the inverse transform is 

Using the Fourier transform, a little thought shows that the energy of a particle in a state 
ip in -L^(A) is 

(cp,hAip) = -^ f dk\ip{k)\''ek, (2.7) 

(27r)'* Ji^^^^u 



with Ek defined in (|1.2| ) and with ip{x) = if 2; ^ A in (|2.5| ). 

Let us consider N orthonormal functions ipi, . . . , ip^, and let E\(ipi, . . . , (p^) be their 
energy. We have 

EAi^Pi,...,ipN) = Y7rTd dkp{k)ek, (2.8) 

with 

N 

p{k) = Y,\'Pjik)f- (2.9) 

3=1 

The function p{k) satisfies the following equations: 

^ p{k) ^ |A|, (2.10a) 

1 



dkpik) = N. (2.106) 
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Indeed, positivity of p is immediate and the last equation is Plancherel's identity. The 
upper bound (2. 10 a) for p{k) can be seen by writing 

p{k) = {f,Pf), (2.11) 

where P is the projector onto {99^}^^, 

N 

P.,y = Y.^j{x)^*{y), (2.12) 

and / is the vector 

/x=e-^'=^XA(x). (2.13) 

Then, since P ^ 1, we have p(A;) ^ ||/||2 = |A|. 
Clearly, we have the lower bound 

Ea,n > inf , ,777^ / dkp{k)ek. (2.14) 



We can use the bathtub principle ( |LL(| , Theorem 1.14) to find the infimum: it is given by 
the function 

I U otherwise, 

where the Fermi level ep is given by the relation , ^ J^, ^^ dfc = A^/|A|. Thus the right 
side of (|2.14D is precisely equal to |A| e(A^/|A|). 



(27r)d Jefe<ep 



3. Lower bound involving the boundary 

In the previous section, we showed that -Ea^at is bounded below by its bulk term. Now 
we strengthen this inequality and prove that -Ba.a^ also includes a term proportional to 
the boundary of A. This can be checked for d = 1 by explicit computation, but higher 
dimensions require more elaborate treatment. 

We start with a lemma that applies when the density n is small enough (or, by the 
symmetry for h\, when it is close to 1). 

Lemma 3.1. //n ^ |5d|/(47r)'^, we have 
Proof. Recall that 



with p{k) = X]i=i l<^(^)P- We want to show that p{k) cannot be too close to pmm{k) in 
( 2.15| ). By completeness of the set of eigenvectors {^j}, we have 

|A| 
pik) = \A\- Yl mk)\'- (3.2) 

j=N+l 



10 
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We now use the Schrodinger equation. We have h\ipj{x) = ejipj{x) for x G A. Let us 
take ^j{x) = for x ^ A; then the following equation holds true for all x G Z : 



■ 22 ^3 {x + e) + Xk'^ {x) 22i '^i {x + e) + 2dipj (x) = ej(pj (a 

e eix+eeA 



(3.3) 



The middle term of the left side is necessary for the equation to hold at sites outside A, 
that are neighbors of sites in A. Taking the Fourier transform, we get 



£Wj{k) + {bk,(fj) = ej(pj{k), 
where b^ is a 'boundary vector', 

bk{x)=X9A{x)e-"''' Y^ e-^^^ 

e:a;+e^A 

Notice that I (9A I ^ ||6fc||2 ^ (2df\dA\ ii \k\oo ^ f. From (O), we have 



(3.4) 
(3.5) 



\ip,ik)\' 



l(^fc,V?j)|^ 
{e,-e,y ^ (4d) 



:^ 



■\{h,Vj)\' 



(3.6) 



The electronic energy in A is given by f p{k)ek- We saw in Section |2| that ^ p{k) ^ 
|A|. By (3^) and (l3.6|) , this can be strengthened to 



const ||P_6fe|p ^ p{k) ^ |A| -const ||P+6fc|p, 



(3.7) 



where P_ (resp. P+) is the projector onto the subspace spanned by (ipi, . . . ,ipj\f) (resp. 
(c^AT+i, . . . , ^\A\))- See Fig. |l| for intuition. We show below that 



\\P+bkf ^ const |aA|, 
and this will straightforwardly lead to the lower bound. 



(3.8) 



I const|lP+6fc| 




const ||P_6fcf 



Figure 1. Illustration of the expression ( |2.g| ) for E\^n; p{k) satisfies more 
stringent estimates than those stated in ( p.lOal ), and this plays an impor- 
tant role in deriving the lower bound. 
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In order to see that the boundary vector has a projection in the subspace of the eigen- 
vectors with large eigenvalues, we first remark that 

{bk,hAbk)= Yl \hix)-bk{y)\^ ^ \\bkf. (3.9) 

{x,y}:\x-y\=l 

We used the fact that each site x of dA has at least one neighbor y outside of A, and we 
obtained an inequality by restricting the sum over such pairs. Let us introduce N' such 
that ctv' ^ ^ and cat'+i > |. We first consider the situation where N ^ N'. Using first 
Cj ^ 4(i and then the previous inequality, we have 

|A| |A| 

M Yl \ibk,^j)\^ > Y. \(bk,^j)\''ej 

j=N'+l j=N'+l 

N' 
> \\bkf-Y\(^''^^^)\'^J ^ h\bkf- (3.10) 

For I /c loo < f , (H) and ( |3lO| ) imply 

2(4d)S 

We can write a lower bound by proceeding as in Section |2|, but using the bound ( p. Ill) 
for p{k), instead of |A|. The bathtub principle then gives 



Pik) < IAI-TT773W- (3.11) 



EA,N-\A\e{n) ^ —— I dfce,(|A| - -1^) - -i- / dke^ ^^^^ 



(3.12) 



where we introduce e'p such that 



"-(^(i^i-w)!,./'- p'^^' 



-F 



Notice that for n ^ \Sci\/{^T^Y , we have e'p < i, so that e^ < e'p implies |fc|oo < f • This 
justifies the use of ( |3.11| ). We bound the first integral of ( 3.12 ) using e^ > £f, and we 
obtain 

EA,«-|AKn) > m.^^l^^Jk(e,-s,). (3.14) 

One can derive a more explicit expression for the lower bound. First, 



/ dA:(eF - efc) ^ ^ep / 



dk. (3.15) 



£f 



Second we use 1 — ^ ^ cos 9 ^1 — ^0^, to get 



^\k\' ^ Ek ^ \k\'. (3.16) 



vr 



One can use the upper bound of ( |3.16| ) to get 

/ dfc ^ \Sd\C^eFf/\ (3.17) 
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Recall that \Sd\ is the volume of the unit sphere in d dimensions. The lower bound of 
( 3.16| ) allows to write 

25^2/d 

.. > ^^. (3.18) 



Then one gets the bound 



24+2d^l+| 

dkisF-Sk) ^ ,/^ . (3.19) 



' ek<e-p 

Hence the boundary correction to E/^^j\f is bounded below by a(n) |9A| with 



Recall that we supposed N ^ A^', where A^' is the index of the largest eigenvalue that 
is smaller than ^- Were it not the case, we can write, with n' = N'/\A\, 

N' N 

^A,iV = J^ej-+ Yl ^J ^ |A|e(n') + a(n')|M|+i|A|(n-n'). (3.21) 

j=l j=N'+l 

We used the previous inequality to bound the first sum, and Cj ^ | for the second sum. 
This is greater than |A|e(n) + a(n)|9A| provided 

e{n') + «("')^ + ki^-n) > e{n) + a(n)^. (3.22) 

A sufficient condition is that -^n — e(n) — a{n) is an increasing function of n. Computing 
the derivative (the derivative of e(n) is £-p{n), that is smaller than {2tt)'^v?''^ /\Sd\^''^ using 
( 3.16| )), and requiring it to be positive leads to the condition 



^ ^ ..,. .. 2^-2' . ....,. - (3.23) 



\Sd\ 
{2(2vr)2 + 55(l + |)}'^/2- 

The right side is greater than |5d|/(47r)'^. D 



It may seem obvious that the extra energy due to the presence of the boundary increases 
as n increases, until it reaches g- But we can provide no proof for this, and hence we need 
a new derivation for the lower bound with higher densities. We proceed in two steps. First 
we give a lemma that works when the boundary has few nearest neighbors; the proof is 
similar to that of the previous lemma. Then we give three lemmas, with more intricate 
demonstrations, and that establish the lower bound for boundaries where at least a density 
of sites have nearest neighbors. We need some notation to characterize the configuration 
around a site x of the boundary. 

Let e, e' be unit vectors in Z ; the notation e\\i means that e is parallel to the z-th 
direction; equivalently, the components of e are given by e^ = iiJifc- We introduce integers 
qx^i and qx,ij', for x € (9A, we set 

qx,i = #{e Wi-.x + e^K] 
q^ij = #{{e, e) : e \\ i, e \\ j, x + e e dA, x + e + e ^ A}. (3.24) 

Notice that ^ qx^i ^ 2 and ^ qx^ij ^ 4. Also, qx^u = #{e \\ i,x + e G dA,x + 2e ^ 
A}, and ^ qx,ii < 2. We also define qx = Yli(lx,i- 
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The following lemma applies to domains where most boundary sites x satisfy qx^ij = 0, 
in which case x has no neighbors that belong to the boundary. Here, qx^ij = means that, 
at X, qx^ij = for all 1 ^ i,j ^ d. 

Lemma 3.2. For all KdZ'^ with 

#{x G M : g,,,, ^ 0} ^ ^|5A|, 
there exists a{n) > such that 

Ea,n ^ |A|e(n)+a(n)|aA|. 
Remark: lim„_>o ctin) = and a(l — n) = a{n) by the symmetry for h^. 

Proof. We can suppose N ^ '-^. The definition of p{k) involves a sum over the first A'" 
eigenvectors (more precisely, of their Fourier transforms). In case of degenerate eigenvalues 
one is free to choose any eigenvectors. For the proof of Lemma 3^ it turns out that the 
possible degeneracy of 2d brings some burden, and it is useful to redefine p{k) by averaging 
over eigenvectors with eigenvalue 2d: 



Ef=l IMk)? + T^^ E,:e,=2. I^.(^)P if A^ > ^V; 



p(fc)=<;;:]v^V;.;;2 , JV^TLv '-^'-^'^ :r.r.\~r. (3-25) 

|A|-2iV ^j-ej = 



here, A" is such that e^ < 2d and e^_,_-^ ^ 2d. The degeneracy of 2d is |A| — 2N (which 
may be zero). Of course, -Ea^at is still given as the integral of p{k) multiplied by e^. 

The goal is to prove that p{k) cannot approach p^i^ in (|2.15| ). Since Ylj=i l^ji^)]"^ — \M^ 
we have 

p{k) ^ |A|-{ Y. \Hk)\' + ^ E \Hk)\'}. (3.26) 

j=\A\-N+l j:ej=2d 

We introduce 

IM 
S{k)= E K^^'V^^OlV- E ■(^^•'V'i)!', (3.27) 

j=\K\-N+l j:ej=2d 

with bk the boundary vector defined in (^.5|). By the inequality (|3.6| ), it is enough to show 
that S{k) is bounded below by a quantity of the order of \dA\. We have 

S{k) = {bk, P+ hk) + Ubk,Po bk), (3.28) 

where P+ is the projector onto the subspace spanned by all ipj with Cj > 2d, and Pq is 
the projector corresponding to the eigenvalue 2d. 

We want to show that S{k) ^ const|5A| for small \k\. This amounts to prove that the 
vector bk cannot lie entirely in the subspace spanned by {(/9j}i ^ j ^tv- 
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Notice that if x G dA and qx,ij = 0, then x has no neighbors in dA. Using the assumption 
of Lemma 3^, as well as |&fc(3;)| ^ 2(i and \dA\ ^ ||&fe|P) we get 

{x,y}:\x-y\=l 

^ 2dY,\bk{x)\''-2d Y, \bk{x)? (3.29) 

xedh x£dA 



The last inequality uses the assumption of Lemma |3.2| , and the fact that |^fc(2^)| is at most 
Id and at least 1. 

Next we consider ||(/iA — 2(i)6fc|p. We have, for x G A, 

[(/iA-2ci)6fc](x) = -^6fc(x + e), (3.30) 

e 

and therefore, if |/c|oo ^ f , 

||(/iA-2d)6fcf ^ ^^|6,(x + e)|2 = ^(2d-<7,)|6fc(x)|2. (3.31) 

a;eA e x 

We write bk = b' + b" , with b"{x) = bk{x) if qx = 2d, otherwise. Notice that b' ± b" . 
Clearly, PqI" = b", and therefore 

Sik) = (6', P+ b') + 1(6', Po b') + i||6"f. (3.32) 

Furthermore, from ( |3.2g| ) and ( |3.31| ), b' satisfies 

ib',{h^-2d)b') ^ -j^Wbkf, {b',{hK-2dfb') ^ ||6'||2. (3.33) 



Because \ej — 2d\ ^ 2d, the last inequality implies 

-Y^\{^^^b')\\e,-2d)+2 ^ \{^^,b')\\e,-2d)=Y,\{v,.b')W^-2d\ ^ 



2d 

j:ej>2d j 



(3.34) 



Y: \{^,,b')\\e,-2d) ^ ^-M!, (3.35) 



With the first inequality in ( |3.33| ), this yields 

Ad 8d 



j:ej>2d 

hence 

E !(»=.. 'V>||^-^. (3.36) 

j:ej>2d 

Back to ( p.32| ), we obtain 

ll/i'l|2 ll/,l|2 ll/,l|2 

5(^) ^ ^^ + 2llHI -J^^ J^- (3-37) 

We can combine this bound with ( 3.26 ) and ( |3.6|) ; we have then for all |/c|oo < f 

p(fc) ^ |A| - ^. (3.38) 
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We introduce ef such that 

1 /, ., |aA| 



(27r)^ 
and we have 



M 



(MY 



dk + 



£fc<ep,|A;|oo<f 



|A| 
{27rY 



dk = N, 



(3.39) 



efe<eF,|fe|oo>-T 



i5;A,^ - |A|e(^) ^ 7T^^(|A| 



(27r)« 



|gA| 
(4d)4 



dfcefe 



+ 



|A| 



ep<£fc<ep,|fc|oo<f 

1 \dA\ 



dksk 



dkek- (3.40) 



We bound the first two integrals using Sk > £f] from the definitions of ep and ep we have 

|9A| f ,,_ |A| 



dfc — 
As a result, we obtain the bound we were looking for, 



dk. 



£F<£k<^F 



i?A,^ - |A|e(^) > 



\dA\ 1 



1^1^ ^ iMY{2TrY 



efe<ep,|fc|oo<^ 



d/c (ep - Sfc) • 



(3.41) 



(3.42) 



D 



We present now another lemma that claims the lower bound for E^^n, and that involves 
a new assumption. We shall see below in Lemmas 3.4 and 3.5 that for all volumes, at least 
one of these lemmas applies. 

Lemma 3.3. Let 6 > and n ^ \Sd\/{^T^Y ■ We assume that 

\\{hK-eN)hk,f ^ <5|9A|, 

for some ko belonging to the Fermi surface, i.e. e^^ = ep where e-p is the Fermi energy for 
density '^ = m- Then we have 

^A,iv - |A|e(n) ^ r/|aA| 

with T] = |5^|5530d+2/(2271d+23^10d+2^130d+9)_ 

The constant ij that appears as a lower bound seems ridiculously small, but we cannot 
do better. Notice that this lower bound is much smaller than the one obtained in Lemma 
3.1 at low density, with n = |S'rf|/(47r)^. We expect however that the lower bound is an 
increasing function of n for ^ n ^ ^, although we cannot prove it. 



Proof. We have 



Ea,n - |A|e(n) 



(27r)^ 



dk 



[— 7r,7r]' 



A_{k)ek-A+{k)ek 



(3.43) 



where 



A+{k) = {\A\-p{k))x[ek<eF] 
A_(/c) = p{k)x[ek > ep]. 
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We used here another convention for the characteristic function, namely X [•] is 1 if • is 
true, and is otherwise. Notice that j dk l\-.{k) = J dk A^{k). Then we both have 

f TTXd f d/c A_ (k) (sk - ep) 
^A,iv - |A|e(n) ^ r^T^ ;,,. , (3.44) 

l(2^i dA:A+(/!:)(eF-efc)• 
And by Holder, this implies 

EA,N-\A\e{n) ^ (^^ J dk[A±{k)]^/^)y [-^ J dk\ek - 8^]--^)' . (3.45) 

One shows in Lemma A.l| (a) that the integral of |efc — efI^^ is bounded by 2. 

Recall that {(pj}i ^j^ \a\ are the eigenvectors of /ia- Let P_ , resp. P+, be the projectors 
onto the first A^ eigenvectors, resp. the last |A| — N eigenvectors. By ( |3.6D , one has 
inequalities 

A-(A;) ^ J—\\P_bkf ifefc>eF. (3.46) 

Let us introduce sets A and A' by 

A = {k : Ek < ep and \k - ko\ < ^25^572 }> 
A' = {k:£k> ep and \k - ko\ < ^25^25/2 }■ (3-47) 

We obtain a lower bound by substituting ( p. 46 ) into ( |3.45| ) , and restricting the integrals 
to A and A' . Namely, 

^A,A^ - |A|e(n) > ^{ji^, J/k\\PM^/')\ 
i^A,^-|A|e(n)^^(^/^^d.||P_5,f/f. 

Let bk = bk/\\bk\\- From the assumption of the lemma, and using ||/ia — e|| ^ 4d and 
Lemma |A.1| (d), we have that for all k ^ AU A' , 

Extracting a factor |i9A|, and using the above inequality, we can write 
^A,^-|A|e(n) ^ \dM ^,,J^^^^,, {l/k\\PM'^^y, 

EA,N-\A\e{n) ^ I^A] ^^3^^^^^^^, (|^^ dfc ||P,b,f /^)'. (3.49) 

Consider k £ A. The assumption of the lemma for kQ, together with the bound for the 
gradient in Lemma |A.1| (e) , implies 

(6fc,(/lA -eAr)^6fc) . 5 



\dA\ 
Therefore 



^ -. (3.50) 



{bk,{hA-eN?h) ^ ^. (3.51) 
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This can be rewritten as 

Y^\{^^,bk)\\e,-eN)^ > ^, (3.52) 

that is, 

Y,\{^j,ik)We] + e%) ^ — + 2eN{h,hAbk). (3.53) 

Hence 

|A| g2 . 

The quantity in the brackets is negative for j ^ N. Observing that e^ ^ siUn)'^ ^ 

1/2 "'"^TT^ (because n ^ |5d|/(47r) and using Lemma AJ (c)), the bracket is bounded by 
2a!+4^2^2_ Therefore, 

{bk,hAbk) ^ eiv + 2'^+W||P+6,.f-^. (3.55) 

On the other hand, for k' £ A' , 

|A| 
(bk'jhAbk') > E Kv'iibfcOPej ^ CAT+i -4d||P_5fc/|p. (3.56) 

j=Af+l 

Since \k - k'\ ^ 6^/2'^'^d'^^/'^, we have from Lemma |A1| (g) and ( |3.4^ ) 

(bk',hAbk')- (bk,hAbk) ^ ^g-^g. (3.57) 



Therefore 



e^+1 - 4d||P_ Vf ^ e^v + 2''+\^d^\\P+hf " ^5^ + ^e^- (3-58) 

Clearly, e^ ^ bat+i; then 

We use now ( p. 49 ). The worst situation happens when ||P+&fc|p is equal to the right 
side of the previous equation. Using Lemma A.l (b) we finally get the lower bound of 
Lemma t3.3|. D 



Now we show that we can use Lemma 3.3 for all A such that Lemma |3.2| does not apply. 

Let a^ = [{2d - £)Qx,i)^^i^a and Qx = ((1 + ^ij)'lx,ij)i ^^j ^^^ More generally, we 
let a denote a vector with entries {2d — e)qx,i, and Q a matrix with entries 2qx^ii in the 
diagonal and qx,ij off the diagonal, that correspond to a possible configuration around x. 
With c = (cos ki)i ^i ^d, we introduce 

F{c] a, Q) = (a, c) + iTr Q - (c, Qc) . (3.60) 

This function appears when establishing a lower bound for ||(/ia — e)&fc|P- 
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Let Q be the set of all matrices Q (for which there exists some compatible configuration); 
we introduce 

Q" = {Q G Q : Qii = 2 and Qi^ + Q^i = 4 for all i / j} (3.61) 

and 

Q' = {Q ^ Q" : Qij ^ 0}. (3.62) 

The reason behind the definition of Q' is that we can provide a lower bound only if 
F{c\ a, Q) is not uniformly zero when k moves along the Fermi surface (i.e. with e^ = ef); 
and we can show that F{c;a,Q) is not uniformly zero only for Q & Q', see Lemma |3.5| 
below. 

For given e^, we define 

fi{eF) = min min max \F{c;a,Q)\. (3.63) 

a,QeQ' ee[0,2d] c:Sk=eF 



We state a lower bound involving //(ep), and check below in Lemma |3.5| that //(ep) is 
strictly positive for ep > 0. 

Lemma 3.4. Let d ^ 2. For all finite A satisfying 



#{x G M : g,,,, ^ 0} ^ ^|9A|, 



we have 



,2 ^ ^(^f) 



^m„p,-e)M^>^^|<>M. 



The factor ^t^i^ is arbitrary here, and has been taken such in order to complement the 
condition of Lemma 



Proof. Let us introduce 

qk{x)=XdA{x) ^ e-^'=^ (3.64) 

e:x+e^A 

By the definition of the discrete Laplacian, 

[(/iA - e)bk] (x) = e-'^'^ [{2d - e)qk{x) - ^ e"^^^ q^ix + e)}. (3.65) 

e 

Let us denote by rfc(x) the quantity inside the brackets above. Clearly, ||(/iA — e)^fc|P = 
||rfc|p. Let TZa, o = Ij • • • , 2*^, represents all combinations of inversions of some coordinates. 
We have the following inequality: 

1 . ^ II 1 - , 2 

a=l a=l 

Indeed, starting from the RHS, we have in essence (with ^ aj ^ 1 and Y2i ^i ~ 1) 



^ ii\2 



^ iY^V^iV^MW) (3.67) 

< [(E».)"'(E' 



1/2 /%r^ „^,|2\l/2 



, fliFil 
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which is the LHS of (|3^ ). 

The RHS of ( |3.66 ) is clearly smaller than max^.^^^^p II (^A — £)^fe|P- One computes now 
Y^a=i ^Tlakix) for X £ dA. First, 

2'^ d 

-^ ^(2d - e)qnAx) = C^d - e) ^ qx,i cos ki. (3.68) 



2'i 

a=l i=l 



Second, 

2d 



4E E ^-""^'^ E ^ 



, , „ , „ i7^aA,•e' 

2^ 

a=l e:x+e£dA e':x+e+e'4A 



d ^ 

'^qx,iiCOs{2ki) - - ^ qx,ij[cos{ki + kj) +cos{ki- kj)] (3.69) 



2 

2=1 «j:«^i 



-Z y ^ Qx,n COS fCj + / ^ Q'xjM / ^ Q'a^iU COS /Cj COS kj . 



We used cos(2A;j) = 2 cos ki — 1, and the bracket in the second line is [•] = 2 cos ki cos kj. 
Gathering ( |3.68| ) and ( p.69| ), we obtain 

1 '' 

jd^'^T^akix) = F{c;a^,Qx). (3.70) 



2'i 

a=l 

One can check that whenever Q^ ^ Q' and differs from 0, there exists a neighbor y that 
belongs to Q' . Then the condition of the lemma implies that 

#{:e GdA:Q,G Q'} ^ ^|5A|. (3.71) 

Furthermore, Q' has less than 5"^ elements since ^ Qij ^ 4; then for any A that 
satisfies the assumption of the lemma there exists Q £ Q' such that 

#{x GdA:Q, = Q} ^ ^^^|9A|. (3.72) 

We get a lower bound for ||(/ia — £)^fc|P by considering only those sites, i.e. 

max \\{hj,-e)bkf ^ ^ max \F{c;a„Q,)\ ^ T^^f^lMI (3.73) 

x&dA:Qx=Q 

uniformly in e G [0, 2d]. D 

There remains to be checked that /u(eF) differs from 0. 

Lemma 3.5. For all sp > 0, we have ^(ep) 7^ 0. 

Proof. We proceed ab absurdo and explore ways where F(c; a, Q) could be uniformly zero. 

The constraint e^ = ep takes a simple form, namely {u,c) = d — gEp- Furthermore, c 

satisfies |c|oo ^ 1; if ep 7^ 0, we can find 6c such that |c+5c|oo ^ 1 and (n, c+6c) = d—^ep 

— in which case 6c must be perpendicular to u. The condition F(c+6c; a, Q) = F{c; a, Q) 
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for all 6c -L u implies that a — 2Qc \\ u. This should also be true when c is replaced with 
c + 6c, hence Q6c \\ u for all 6c _L u. Now take {6c)i = Sn — 6ji. We have 

{Qdc)j = -Qjj + Qj„ (3.74) 

and these two components must be equal, since Q6c is parallel to u. Hence Qa + Qjj = 

In this case F(c; a, Q) takes the form 

d d d 

F{c; a, Q) = (2d - e) ^ qx,iCi - {2d - e-p) ^ fe.iiQ + ^ Qx^u- (3.76) 

Since Q £ Q' we have g^^y ^ 0; if (u, c) = 0, one can take c = 0, and F(c; a, Q) is 
strictly positive, so we can suppose {u, c) / 0. 

Let s = J2i Qx,ii/iu, c), and w the vector with components 

Vi = {2d - e)qx,i - {2d - eF)qx,ii + s. (3.77) 

Then F{c; a, Q) = {v, c). If we require this to be zero for c || u, then we need v -L u. But 
we also require {v, c + 6c) = {v, c) for all 6c _L n, hence v \\ u. So v must be zero, i.e. 

{2d - e)q^^i - {2d - eF)q^,ii + s = (3.78) 

for all 1 ^ i ^ d. 

We also have qx,i + Qx,ii ^ 2, and q^^u cannot be always equal to 2. If s / 0, one checks 
that necessarily q^^u = 1, which is impossible because Q G Q' . Hence F{c; a, Q) cannot be 
uniformly zero when moving along the Fermi surface. D 

4. Finite U 

We consider now the Falicov-Kimball model with finite repulsion U, and establish a 
lower bound for the ground state energy of N electrons in a configuration specified by A. 
More precisely, we show that when decreasing the repulsion U, one does not lower the 
energy more than const • \dA\/U. 

For any A, the spectrum of h^^ is included in [0, 4(i] U [U,U + Ad]. When U > Ad, 
eigenstates with energy in [0, 4d] show exponential decay outside of A; and eigenstates 
with energy in [U, U + Ad] show exponential decay inside A. Hence A and $7 \ A are 
essentially decorrelated, and the situation is close to that with U = oo. 

The following proposition compares the energies of electrons with finite and infinite U. 
It is useful to introduce r]{U), 

/.n ( 2d \^4^\{U-2dfy {U-2df'^ ^ .^. 

^(^) = iu^) g [wr^)\ = {U{u-L)Y - '■ ^'-'^ 

Notice that lim(7^oo U'^'>]{U) = 4(i^, as it easily comes out from the middle expression. 
Proposition 4.1. If U > Ad, we have 

eIn > E^^N - j{U)\dA\, 
with 

l{U) = jf^ + d2'+'v{U). 
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Proof. First, we remark that eigenvectors of h^ j^ with eigenvalue smaher than 4d have 
exponential decay outside of A. Indeed, for x ^ A the Schrodinger equation can be written 

If ej ^ 4(i, we have 



U + 2d- ej 



Li^.(-)i ^ jjr^^2 ■ (4-3) 

Using this inequality, we can proceed by induction on the distance between x and A. 
The induction hypothesis is that the following holds true 

El-i(^)l^ , (^)=" (4.4) 

J = l 

for any x such that dist (x. A) ^ n. As a result, we have 

V^l ^ m2 / / 2d N2dist(x,A) 

Let us introduce 

AT 

Pxy = ^Xa{x) V*{x) ifjiy) XK{y). (4.6) 

j=i 

We show that E^^ is bounded below by Tr ph\, up to a contribution no greater than 
const \dA\/U. Recall that /ia is the Hamiltonian with infinite repulsions. If Pa is the 
projector onto the domain A, let ipj = Pa^j- 

N 

i=l {x,y}:\x—y\=l x^A 

N 

^ E( E l'^i(^) -'^i(y)l^+ E l'^i(^)-'^i(y)l^ - E l'^i(^)l^ 

i=l {x,y}:|x-i/|=l {z,j/}(Z:A zeA,j/^A 

|x-y|=l |a;-|/|=l 

w 
^ Trp/iA-2^ ^ 1^,(^)1 I (^,-(y) I . (4.7) 

i=i xeA,j/^A 

|z-y|=l 

By the Schwarz inequality, the last term is smaller than 

N , /„ ^ 1/2 0^2 

K E Ei^.(-)ir ( E Y.\^m') ^iT^i^Ai- (4-8) 

x&K,y^Aj=l x&A,y^Aj=l 

\x-y\=l \x-y\=l 

We used ( [4. 5] ) with dist (a;. A) being respectively and 1, in order to control the quantities 
in both brackets. 
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Recall that ej denotes the j-th eigenvalue of the Hamiltonian /ia; that is, with infinite 
repulsions. Let us introduce the projector Pj onto the corresponding eigenstate. Then 

|A| |A| 

Tr /5/iA = 2_^ ^j '^^ P^j = z^ ^i ""■?' ("^-S) 

i=i i=i 

where the rij satisfy ^ rij ^ 1 , and ^ • rij = Tr p. By the bathtub principle PjL| , we 
obtain the lower bound 

Trp 

TVp/iA ^ ^ej. (4.10) 

i=i 
There remains to show that Trp is close to A^. We have 

iV-Trp=5]^|v.,(x)p ^ ^#{x:dist(x,A)=n}(^^) 



x^K j=l n=l 

n=l 

We bounded #{•} ^ 2'^("+f7^)|9A|. Since e^ < 4(i, we obtain the proposition. D 



5. Upper bound 

We establish now an upper bound for the sum of the first A^ eigenvalues in a finite 

A,N 



domain A, for the case of infinite repulsion. The bound carries over to finite U, since E 



is increasing in U. 

The strategy is to average h\ over a huge box. The 'strength' of the averaged Hamil- 
tonian depends on the number of bonds in A, which is roughly 2d\A\ — \dA\. The averaged 
Hamiltonian is, up to a factor, the hopping matrix in the huge box, and its ground state 
energy is easy to compute in the thermodynamic limit. This can be compared to £'A,Ar by 
concavity of the sum of lowest eigenvalues of self-adjoint operators. The result is 

Proposition 5.1. The sum of the first N eigenvalues of the Laplace operator in a domain 
A with Dirichlet boundary conditions, satisfies the upper bound 

Ea,n ^ |A|e(n)-F|aA|(2dn-e(n)). 

Proof Let L be a multiple of |A|, and A''^ be such that Nl/L'^ = A^/|A|. We consider a 
box {1, . . . , L}'^. We introduce e = \{eN + cat+i). Let TZa, a = 1, . . . , L'^dl^ represent a 
translation possibly followed by an axis permutation. We define the averaged Hamiltonian 

L'^d\ 

hL,A = j^^^yhTz^A-ilTZaA)- (5-1) 

a=l 

Then 

L'^dl 

■ a=l 

= SN{hA-etA). (5.2) 
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Indeed, all summands in the above equation are equal, and the Hamiltonian /ia — sIa has 
no more than N negative eigenvalues, and at least L'^ — |A| zero eigenvalues. The RHS is 
equal to -Ea^at — A^e. 

Let Ki be the number of sites in A that have i neighbors in A. We have | A| = X^j^g -^« 
and \dA\ = Yli=o ^i't and the number of bonds in A is 2 X^j=o^-^i- Then the averaged 
Hamiltonian is 

{hL,A)^y = -j^S\^_y\=i + {2d - e)j^6^y, (5.3) 

with 

2d 



2d 

j=0 



^2d 2d 



Let K = ^^^0 ^Ki; then t = \A\ - K and K ^ \dA\. One easily checks that 



Notice that all operators commute. In ( |5.2D , the terms involving e cancel, since 5'Ar^(]lM^ ^^id) 
Nl, and Nl^-^ = N. Now, as L ^ cxd, 

;^'S'Afz.(/i{i,...,L}d) -^ e(n). (5.6) 

Therefore (^^ implies 

\A\e{n) + K{2dn-e{n)) ^ Ea,n- (5.7) 

D 

6. Positive electronic temperature 

This section considers the electronic free energy at positive temperature, for a fixed 
configuration of classical particles. We will see that the inequalities satisfied by the sums 
over lowest eigenvalues have an extension to free energies. 



6.1. Lovi^er bound for U = oo. We start with U —>■ oo. Let Fa{P, /i) = — 4 log Tr e f^^^ 
where the trace is taken in the Fock space of antisymmetric wave functions on A, and 

tU=o 

'n,A 



H\ = HR . °° is the second quantized form of the one-particle Hamiltonian h\ defined by 



Proposition 6.1. For all finite A, one has the lower bound 

F^{P,^^)-\A\f{P,^^) ^ a{(3, fi)\dA\, 

where a{l3,fi) > satisfies lim/3^00 '5(/3, /u) > if < fi < Ad. 

Proof. The fermionic free energy F\{P,^) can be expressed in terms of the eigenvalues of 
/lA, 

1 1^1 
i^A(/?,/x) = --.5]log(l+ e-^(^^-'^)). (6.1) 
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In order to compare this with the corresponding infinite- volume expression ( |1.9[ ), we par- 
tition the BriUouin zone [— 7r,7r] according to the level sets of the function Sk] more 
precisely, we define measures Hj, 1 ^ j ^ |A|, by 

d/^.(^) = (1^ ^ [^f(^) < e{k) < eF{^^)]dk. (6.2) 

Notice that J dfj,j{k) = 1 and -m- X^Ui d/ij(fc) = ,i^^,^ . Next we introduce e*, that are 
equal to Sk averaged over fij-. 



^i 



/ dtij{k)ek. (6.3) 



The ground state energy ( p..3p of a density -/V/| A| of electrons in Z can then be written as 

1 ^ 
e(iV/|A|) = -j;e*. (6.4) 

' ' i=i 

From the lower bound without a boundary term, we have 

N N 

^e,>^e*, (6.5) 

for all A^ < |A|, and equality when A^ = |A|. 

Actually, inequality (|6.5| ) can be strengthened by introducing a term depending on the 
boundary of A. In Theorem |1.1| , a{n) can be taken to be increasing in n for n ^ g- Also, 
a(l — n) = a{n). Therefore there exists a function a{e), with a(e) > for < e < 2d, 

a{Ad — e) = —a{e), and 

"("-) = TTTZw / dka{ek). (6.6) 

'£fc<eF 



(27r)'^ 



:;. = |d^,(A;)(e, + ^a(efe)); (6.7) 



Next we define 

e 
then the following is stronger than ( |6.5D and holds true, 

N N 

Y^e, ^ J^e;.. (6.8) 

i=i i=i 

With a(e) chosen appropriately both sequences (cj) and (e') are increasing, and the in- 
equality above is an equality when A^ = |A|. The sequence (e^) is said to 'majorize' (e'). 
We can apply an inequality due to Hardy, Littlewood and Polya (and independently found 
by Karamata); see [Mit| page 164. For any concave function g, we have 

|A| |A| 

Y,9ie,) > j;<7(e;). (6.9) 

i=i i=i 

(Conversely, if (^]^) holds for all concave g, then (cj) majorizes (e').) We use this inequal- 
ity with 

5(e) = -ilog(l+e-^(^-^)), (6.10) 
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which is concave. We get 

Fa(/3,/z) ^ ^9(e;) ^ -^ / dk9{e,+ \^a{ek)), (6.11) 

where the last step is Jensen's inequahty. Then 

^Fa(/3,/x) -/(/?,/.) > ■^ldk{g{ek+^-^a{ek))-g{ek)}. (6.12) 

In the limit /? — > oo, we have 

.(.j-lr" ''^1" (6.13) 

10 it e ^ //. 

As a result, for aW < fj, < Ad we get a lower bound for large P that is uniform in the 
limit /3 ^ c». 

One also gets a lower bound by using concavity of g, that holds for all temperatures, 
but that is not uniform in /3: 

\9M , , 

1 If f'^k + -rxra{ek) 



^^-F^{P,^l)-f{f3,^l) ^ j^J dkj^ deg'ie) 



L_ fdk{\^a{e,)g'{ek)-0{a{e,f)} 



{2-kY 
\dh\ 1 



dka 



{ek)[g'{ek) - g{M - Sk) - 0(a(efe))}. 



|A| {2^YJe,<2d 

(6.14) 

The integrand in the last line is strictly positive if a{ek) is small enough, and chosen to 
vanish appropriately as e^ -^ 2d. D 

6.2. Louver bound ■with finite U. We extend now the results of the previous section to 
the case of finite repulsion U . As we noted in Section ^, when U > Ad all eigenstates have 
exponential decay, either in A or in Q \ A. We show that the total free energy in is equal 
to a term involving il \ A only, plus a term involving A only, up to a correction of order 

\dK\/U. 

Proposition 6.2. 

FnAf^.tA ^ FAiP,f^) + Fn\A{f3,fi-U)-j{U)\dA\ 
with 

7([/) = {2dU + Ad + 8d^)2'^r]{U) + ^^"^^ 



U -2d' 
and rj{U) is defined in (|]l|). Notice that limjy^oo U^{U) = 16d^ + 8(i^2'^. 

Proof. Let us introduce 

_ J(/jj(x) if 1 ^ J ^ |A| and x G A, or if |A| < j ^ \9.\ and x ^ A 

I otherwise. 
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We assume A^ > |A| (otherwise, replace |A| by N in the next expressions, and ignore the 
sums whose initial number is greater than the final one). Then 

N |A| 

3=1 j=l {x,y} {x,y}glA 

\x-y\=l \x-y\=l 

|2^ 



Y^ \^,{X)\' + UY,\V,{ 



i[X 



x£A,y^A 
\x-y\=l 



xiA 



N 



+ E ( E l'^i(^)-'^i(y)l^+ E l<^i(^) -</'i(y)l' 



i=|A|+l {x,y} 
\x-y\=l 



{x,y}^A'' 
\x-y\=l 



Y^ \^,{x)\' + uY\v,{x)\- 



x^A,y£A 
\x-y\=l 



x4A 



We proceed as in Section |^ and define 



i=i 



N 



j = |A|+l 



Then 



N 



N 



i=i 



j=l xeA,y<^A 
\x-y\=l 



The inequality (4.10) is still valid, for both p and p' . Hence 

N Trp |A|+1+Trp' AT 

E^f ^E^i+ E i^j + u)-2Y E \^jix)\\^3iy)\- 



i=i 



i=i 



j=|A|+l 



i=l xGA,yi^A 
\x-y\=l 



(6.16) 



(6.17) 



Y ef > Tr phA + Tr p' [hn\A + f/] - 2 J] ^ \^j{x)\\^j{y)\. (6.18) 



(6.19) 



(6.20) 



Here, ej, |A| < j ^ |ri| are the eigenvalues of the operator /if2\A- 
We define 

' ^dZx^^\Vjix)\^ + '^ExeA,ytA\V3ix)\\'Pj{y)\ ifl ^ i ^ |A| 

r X — 1j\ — ]_ 

^■" |(4d + t/)E.6Alv5i(^)P + 2E.6A,j;^Al95,(x)||<^i(y)| if|A|<i ^ \n\ 

\x~y\=l 

Then ( |6.19 ) takes the simpler form 

N |A| N 

E^f > E(^^-^^)+ E i^J+U-5,). (6.21) 

i=i i=i i=|A|+i 

The sequence in the RHS is not necessarily increasing, but one gets a lower bound by 
rearranging the terms. Hence one can apply Hardy, Littlewood, Polya inequality. Indeed, 
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it also works when the total sum over elements of the sequences are not equal, provided 
the concave function is increasing — which is the case with g{e). One obtains 

\n\ |A| |n| 

E5(ef) ^ E5(^J-'^i)+ E 9{ej+U-5j). (6.22) 

i=i i=i i=|A|+i 

We use now g{e — 5) ^ g{e) — 5, and we find 

|n| 

FE,A(f^^ ^) > ^a(/?, /x) + Fn\j,iP, f,-U)-Y,^J- (6-23) 

i=i 

The remaining effort consists in estimating the sum of 5j, using exponential decay of 
eigenfunctions ipj either in A or in 17 \ A. Retracing ( 4.11| ) and (4.8), we get 

'^' 8r/2 

^6j ^ Ad2''rj{U)\dA\ + j^-^\dA\, 

'""' 8./2 

5] 6, ^ (U + Ad)2'r^{U)\d{n\A)\ + j^-^\dA\. (6.24) 

i=|A|+i 

Notice that the last term can be written with \dA\ instead of \d{0, \ A)|, as can be seen 
from (U). We use d{n \ A)| ^ 2d\dA\, and we finally obtain 

'^' (4d)^ 

^6j ^ {2dU + 4d + 8d^)2'^r]{U)\dA\ + -^-L^\dA\. (6.25) 

3=1 

u 

6.3. Upper bound. We turn to the upper bound for the electronic free energy. We first 
notice that the free energy is raised when one decorrelates the domain occupied by the 
classical particles, from the empty domain. The following proposition applies to all finite 
subsets of Z , and it also applies when O is a finite d-dimensional torus. 

Proposition 6.3. We have the upper hounds 
. Fs^A(/3,/i) ^ FA(/3,/x)+F^\A(/3,/i-[/). 



Fa(/?,/.) ^ |A|/(/3,/i) + -^(^|A|^+2d|M|). 



, d-l 



Notice that the isoperimetric inequality implies that for all finite A C Z , |A| d ^ 
\dA\. This does not hold, however, when A is e.g. a box with periodic boundary conditions. 

Proof. The Peierls inequality allows us to write 

j 

for any set of orthonormal functions {V'j} (in the Fock space of antisymmetric wave func- 
tions on Q). We can choose the ij)j to be eigenfunctions of H{^ and Hq\j^ — decorrelating 
A and fi \ A. In fi \ A, the free electrons experience a uniform potential U; the energy 
levels are given by the spectrum of /ij^x a plus U. This only shifts the chemical potential, 
so that we obtain the first claim of the proposition. 
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Now we estimate F/^{(3,ij,). Let us introduce 

e- = (1 - If )e* + 2(i^ = I d/.,(fc) [sk + (2d - e,)^-^] ; (6.27) 

then Cj ^ (^j+ij Ylj=i^j — ^7=1^^! ^'^d ^^6 upper bound for the ground state energy 
can be cast in the form 

N N 

Y,e, ^ Y.^r (6-28) 

i=i j=i 

This allows us to summon again the Hardy, Littlewood, Polya inequality, and we get 

|A| 
Fa(/3,/z) < Y^aU df,,{k)[ek + {2d-Ek)^-^]y (6.29) 

i=i 
The derivative of g{e) satisfies 

< g'ie) ^ ^^^ (6.30) 

(recall that e ^ 0). Since the measure fij is concentrated on those k where e^ lies between 
eF(Tjr) and £F(m)5 we can bound ( |6.29| ) by 

^a(/5,/^) ^ (|^^_^_^^/^5(e.) + 3^^(E[eE(A)-eE(^)]+2d|M|). 

(6.31) 

We need a bound for ^Fijjn) ~ ^f(t^); since Vsk = 2{smki, . . . ,smkd), we have 
||Vefe|| ^ 2\/d. Let us take k such that Sk = £F(Tjr), and 6k \\ k such that efc+<5A: = 
£f(^). Then 

£F(^)-eF(^) ^ 2Vd\\6k\\. (6.32) 

If (J/cmin is choscu SO as to minimize the norm of such 6k, we have 
11/" 1 

|A| (27r)^ 7,p(^)<,,<,p(^) (27r)°' 

Combining this inequality with ( 6.32 ), we get 



This leads to the upper bound of Proposition |6.3| . D 

6.4. Proofs of the corollaries. 



Proof of Corollary \1.3{ . Let e~ = —d + h and e^ = —d — /i be the energies per site of the 
all — and all + Ising configurations. A configuration can be specified by the set A of — 
spins. Let B{K) be the set of bonds connecting A and fi \ A. Notice that ^|i3(A)| ^ 
|(9A| ^ |;B(A)|. The partition function of the Ising model can be written as 

Zi,f, = Y, e-^[l^l'^~+l^\^l^^] ^-msm . (6.35) 

AcO 



SEGREGATION IN THE FALICOV-KIMBALL MODEL 29 

Now the upper bound for F^{^{(3,ij) implies that the partition function of the FaUcov- 
Kimball model is bounded below by 

Z^ ^ y^ g-/3[|A|/(/3,M)+|n\A|/(/3,M-(7)] ^-l3CaJdh\ ^-PC'^^^M^ _ ^g^3g^ 

The last factor vanishes in the thermodynamic limit. One then makes the connection with 
Ising by multiplying Zq, by 

exp{/3|J^|(fC,,^ + i[/(/3,/i) + /(/3,/i -[/)])}, 
and by choosing the temperature to be \Cd,^(i, and the magnetic field to be 

/i = 7^[/(A/^)-/(A/^-^)] (6.37) 

(the magnetic field is negative). The other bound is similar, simply replace C^,^ by 
a/2d. ' D 



Proof of Corollary L4. Because ft was assigned periodic boundary conditions, we have 

{Swa:,Wy)n = T7y^\/_^^W:, + ^,Wy + ^)Q- (6.38) 

It is not hard to check that for any configuration w specified by A C 0, one has 

5^<5..+.,«,,+. > \n\-\dA\\x-y\,. (6.39) 

zen 



Then 

\n 



{Sn..,n.y)n > l-k-y|i(W>n- (6-40) 



We need a bound for the last term. The fact is that typical configurations of classical 
particles cannot have too much boundary: K^ is smaller than r = „-?§ k . Indeed, 



Eaco X [\dA\ >r\n\] e-^^KAM 



Eacc e 



-pfKA^^p-) 



2M p-/3|A|/(/3,/.)-/3(|Q|-|A|)/{/3,M-C/) p-(21og2)|n| 

^ -^ -r, srr (6-41) 

e-/3|A|/(/3,/.)-/3(|c|-|A|)/(/3,M-c/) Q-PCaA^,\n\)'hr- ^-HC'^^^M-^ 

Therefore 

The last term vanishes in the limit il ^ Z^, and the term involving jx — y|i vanishes when 

(3 ^ oo. D 
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7. Conclusion 

Our analysis of the Falicov-Kimball model away from half-filling allows some extrapo- 
lations. We expect segregation to survive at small temperature, when both the classical 
particles and the electrons are described by the grand-canonical ensemble, at inverse tem- 
perature P and with chemical potentials Hc and /ie- Segregation is a manifestation of 
coexistence between a phase with many classical particles and few electrons, and a phase 
with many electrons and few classical particles. It is therefore natural to conjecture the 
following, for d ^ 2: 

A first order phase transition occurs at low temperature, when varying the chem- 
ical potentials. 

The transition from the chessboard state at half-filling for the itinerant and heavy 
electrons (and large U) to the segregated state is still not clear. A heuristic analysis 
suggests that these states could coexist, hence there could be another first-order phase 
transition. Alternate possibilities include mixtures between other periodic phases and the 
empty or full lattice before the segregation sets in. 

One interest of the Falicov-Kimball model is its possible relevance in understanding 
the Hubbard model, a notoriously difficult task. See e.g. [|Lieb2 | and | Tas2 ] for reviews 



of rigorous results on the Hubbard model. The relationship between the Falicov-Kimball 
model and the Hubbard model is like the one between the Ising and Heisenberg models 
for magnetism. The former does not possess the continuous symmetry of the latter, and 
therefore the approximation is a crude one. Still, the two models share many similarities; 
for instance, the Falicov-Kimball model displays long-range order of the chessboard type 
at half- filling and at low temperature ||KL[|, and the ground state of the Hubbard model 



is a spin singlet [Lieb 



Ferromagnetism in the Hubbard model depends on the dimension, on the filling, and 
on the geometry: it has been shown to occur on special lattices such as 'line-graphs' 
| MT| , Tas| , |Mig| , Tas2]. Does ferromagnetism take place in the Hubbard model on Z^, for 



large repulsions and away from half-filling? 

Returning to Falicov-Kimball, let us walk on the road that leads to Hubbard. We 
consider the asymmetric Hubbard model that describes spin 2 electrons with hoppings 
depending on the spins (this interpretation is more convenient than physical). Its Hamil- 
tonian is 

Ht = - Yl 4tST-^ Y1 4iSi + f^2^™^T^xi- (7.1) 

x,y-\x-y\=l x,y:\x-y\=l x 

Notice that Hq is the Falicov-Kimball model, while Hi is the usual Hubbard model. Al- 
though we did not prove it, it is rather clear that segregation still takes place for very 
small t. Furthermore, the density of the phase with classical particles, in the ground 
state, should still be exactly 1 — indeed, the electrons exert a sort of 'pressure' that packs 
the classical particles together, and the tendency of the latter to delocalize is not strong 
enough to overcome this pressure. This is summarized in the following conjecture: 

For t ^ to, segregation occurs in the ground state, at large U and away from 
half-filling, in the form of a coexistence between a phase of classical particles 
with density 1, and a phase of electrons with smaller density. 

This should also hold at positive temperature, although the density of the phase of the 
classical particles will be reduced, due to the presence of some holes. 
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If we increase t, assuming that segregation remains, we should reach a critical value tc < 
1 where the region of classical particles starts to grow. The density of the phase of particles 
with smaller hoppings is now strictly less than 1. A major question is whether segregation 
survives all the way to the point where t reaches 1 — this would imply the existence of a 
ferromagnetic phase in the Hubbard model. We note, however, that while it is conceivable 
that there is a segregated (i.e., ferromagnetic) ground state at t = 1, it cannot be true that 
every ground state (for equal number of up and down spins) is segregated. This follows 
from the SU(2) symmetry. If ^ is a saturated ferromagnetic ground state with 2N up 
electrons, we can construct $ = (S"-) \I', which is also a ground state, with N up and N 
down electrons. However, $ has the up and down electrons inextricably mixed, which is 
the opposite of a segregated state. Indeed, the SU(2) symmetry is restored precisely at 
t=l, and the ground states have at least the degeneracy due to this symmetry. 

The Hubbard model is a rich and complicated model that poses difficult challenges. 
The Falicov-Kimball model can be of some help, for instance in checking scenarios that 
should apply to both models. This discussion of ferromagnetism illustrates however that 
the links between them are subtle. 

Appendix A. 

We derive in the sequel various expressions that are too intricate to appear in the main 
body of this paper. 

Lemma A.l. 

(a) ^Jdk\ek-eF\-'/^<2. 

(b) Assume that c? ^ ]t n'Jd i^^^'^i then for all k such that Sk = £f, 

dk' X[\k'-k\<a] ^ |5d|(^)^. 

(c) e(n) ^ 12(^)V+l/|5,|i. 

(d) IIV^II ^ M^l\ 

(e) ||v H%^f^ll' || ^ 29d"/2. 

(f) livi^^^f^ll ^ 32d7/2. 

(g) Assume that pkf' l\dk\ ^ ry. Then if ri ^ 1, ||V%J^|| < rfH^d^^l'^. 



Proof of Lemma A.l (a). Setting Y = 2d — ep — 2 X^j=2 ^°^ ^*' ^^^ making the change of 

1 



variables ^ = cos ki , one gets 

[ dk\ek-£Y\~^/^ = 2 [ dk2...dkd r dki 

J Jf-7r,7rld-l Jo 



-7r,7r ' 



y-2cos/i:i|V4 

2 / dk2...dkd I dC 

^l — ^^ \r — z^ !-"/■* 

3/4 



1 7-1 \/r^\Y-2c\y'' 

1/3 
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The integral over ^ can be split into one running from to 2) and one running from ^ to 
1. For the first part we bound ^, j*; ,3^4 ^ 2'^/^ -7=, while the bound for the second part 
can be chosen to be V^tt^jtsji- Everything can now be computed explicitly, and we find 
231/1232/3 (2^)^^-! < 2(27r)f D 

Proof of Lemma A. j| (b). Let us introduce a map 7(^) such that 1 — ^7^(0 = cos^; pre- 
cisely, 

|-V2(l-cosO if^G[-7r,0]. 



The condition e^ < ep becomes J2i=i l7(^«)P < ^F- The derivative of 7 is 

d7 |sin^| 



(A.2) 



de v^2(l-cose) 
We check now that \'y{C')-^{S,)\ > !■?' -^P/47r. Let us assume that 7(^') > 7(^). Then 

7(C') - 7(6 = r dA 'f"^' ^ i f dA| sin A| (A.3) 

J( V 2(1 ~ cos A) ^ J^ 



^ - TdAIAI ^ \e-e/^^- (A.4) 

^ Js 



Then we can write 



I dk' x[\k'-k\<a] ^ f dk' x[\k'i-ki\<f2^ 

„ d 

^ / d7i . . . d7^ X [Y^ Wf < ep] X [h- - 7i\ < 



i-irdi 



One gets a lower bound by replacing the last characteristic function by the condition 

J2i=i Wi ~ 7«P < ^- Recall that ep ^ , p/d ?^-^ i the assumption of the lemma implies 

2 
that y^ep > ^^; as a consequence, a lower bound is the volume of the sphere of radius 

oL n 

8wd- '-' 



Proof of Lemm,a A.l (c). By (3.16), 



''"'*^t<./'='"^^^WT^^*"- <"' 



The lower bound then follows from 






Proof of Lemm,a A.]\ (d)-(g). Since 



xe9A e:x+e4A 
e'-.x+e'^A 
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we have 



V^ 



|9A|' 



\dA\ ^ ^ 

' ' xeaA e,e' 



iiej-e'j)e'^^^-^'^ 



(A.8) 



This is less than 2(2(i)^, and we obtain the bound (d). 
We consider now ||(/ia — £)&fc|P- 



e"'-[{hA-e)bk]{x) = {2d-e) Yl «' 



-ifce' 



E 



-ik{e+e') 



e':x+e'^A 



e:x+e(^dA 
e':x+e+e'^A 



Y^ g-ifc(e+e') ((2d _ e) X [|e| = 0] - X [|e| = 1] ) . (A.9) 



e:a;+ee9A,|e|=0,l 
e':x+e+e'^A 



In the last line, e is allowed to be 0. Let ^(e) = {{2d — e)x [|e| = O] — X [|e| = l])- Then 



[(/lA - e)bk\{x)\ 



E 



E 



Jk(e+e'-e"-e"') 



i{e)i{e"). (A.IO) 



e:a;+ee9A,|e|=0,l e":a;+e"G9A,|e"|=0,l 
e':x+e+e'^K e"':x+e" +e"'<^A 



One computes now the j'-th component of the gradient; it involves a term e,- + e'- — e',' — e'" 
that is smaller than 4; there are sums over e',e"', with less than {2d)'^ terms; the sum 
Se 1^(^)1 i^ bounded by Ad; finally, the number of sites where (/ia — e)hk differs from is 
bounded by 2(i|9A|. As a result, the j-th component of the gradient is bounded by i(4d)^, 
and we obtain (e). 

We estimate now the gradient of (6^, h\bk)- One easily checks that 

jk{e-e'-e") 



(A.ll) 



ibk,hAbk) = \\hf-Yl E E E e- 

xedA evx+e^A e':x+e'&dA e" -.x+e' +e" ^A 

We can use the bound (d) for the gradient of ||6fc|p. The gradient of the last term is less 
than 3(2(i)^|(9A|, so we can write 

,(5fc,/iA5fc) 



V- 



\dK\ 
Finally, one easily checks that 



^ 8d5/2 + 24d7/2 ^ 32d7/2. 



(A.12) 



V 



{hk^h^bk) 2 ^ ^f \dK\ ^2 ^{bk^hxbk) 2 , ^( {bk^h^bkW r \dK\ 



\b,P 



< 2 



I6.||2 



V- 



|M| 



+ 2 



\dK\ J \\M\ 



V- 



\b,r 2 



|9A| 



Using (d) and (f), as well as {bk,h[J)k)/\dK\ ^ 2{2df' , one gets (g) 



(A.13) 
D 
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